clear
global lambda1;
options = optimoptions('fsolve','Display','iter','MaxFunEvals',2400); % Option to display output

i=1;
for lambda1=0.05:0.02:0.88    
    [x,fval] = fsolve(@myfun1,[1;3;0.6;1;1.2;0.3],options); % Call solver
    lambda2(i)=x(3);
    i=i+1;
end

t=0.05:0.02:0.88;
plot(t,lambda2),hold on
plot(lambda2,t),hold off